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Abstract. In this paper we study the equidistribution of expanding 
horospheres in infinite volume geometrically finite rank one locally sym- 
metric manifolds and apply it to the orbital counting problem in Apol- 
lonian sphere packing. 

1. Introduction 

In this paper, we study the equidistribution of expanding horospheres in 
infinite volume geometrically finite locally symmetric rank one manifolds 
with respect to Burger-Roblin measure. As an application we apply it to 
the orbital counting of geometrically finite groups. 

A priori it is not clear how to count the growth of the number of or- 
bit points in F""*"^ under the infinite co- volume group T C 0^{n, 1) where 
Ow{n,l) = 50(71, l),S'C/(n,l), 5^(71,1) depending on F = M, CM. KQfN) = 
for a signature (n, 1) quadratic form Qf, we give a quantitative estimate 
of the asymptotic growth 

#{v G ?;or| \\v\\ < K} 

for any norm 1 1 • 1 1 on F"+^ with the control of the error term. Controlling the 
error term is crucial to our application to the counting of prime curvature 
spheres in Apollonian sphere packing. This orbital counting theorem follows 
from the equidistribution of expanding closed horospheres in geometrically 
finite groups. First we show that for any tp S C^{r\G)^ , the average of 
this function on a horosphere of height y can be explicitly estimated in terms 
of L^-product of tp and ^q, and y^~^ where po is a unique eigenfunction of 
the Laplace operator with L^-norm 1 of eigenvalue 6{D — S), and 6 is the 
critical exponent of F, D is the Hausdorff dimension of the ideal boundary 
of the associated symmetric space. See Theorem 19.21 

The techniques involve the unitary representation theory, measure theo- 
retic approach in algebraic Lie groups, Patterson-Sullivan measure on limit 
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sets and some geometrical insights in rank one space. We carry out the 
computation in exphcit coordinates, so-called horospherical coordinates in 
rank one spaces. See sections [2] and El We outline necessary backgrounds as 
the proof evolves in coming sections. 

Let X be a real, complex or quaternionic hyperbolic space with curvature 
between —4 and —1, and F C Iso{X) = G a geometrically finite group whose 
critical exponent is 5 > D/2 where D is a Hausdorff dimension of dX. Let 
G = KAN be a fixed Iwasawa decomposition introduced in section [2j In 
other words, is a (generalized) Heisenberg group, j4 is a one-parameter 
subgroup stabilizing a chosen geodesic (0, 0, y) in horospherical coordinates 
(see section [21 [3]), and K is a maximal compact subgroup stabilizing (0, 0, 1). 

The main theorem is: 

Theorem 1.1. Let T be a geometrically finite discrete subgroup in SO{n, 1), 
SU{n, 1), Sp{n, 1) with the critical exponent (5 > y where D is the Hausdorff 
dimension of the boundary of the associated symmetric spaces of the sectional 
curvature between —4 and —1. Suppose vq is in the light cone such that vqV 
is discrete, and the stabilizer of vq in g^^Tgo is in NM. Then for any norm 
II • II on F"'i 

#{v G ^;or : \\v\\ < T} ~ c^,S-'T^ [ \\vo{go^kgo)\\-^dk. 

Jk 

If \ \ ■ \ \ is g^^K go-invariant, then 

G ^;or : \\v\\ < T} = c^„<^-1t^| |^;o| r^(l + 0{T-^')). 
Here 5' depends only on the spectral gap. 

This theorem is a generalization of |23j to a general rank one case. 

Notation: f{x) = 0{x) implies that liuix^oo < oo. 

Outline of a proof: One introduces a continuous counting function 
Ft{9) = T.'ye{rnNM)\rXBT{voig), where Bt = {v e vqG : \\v\\ < T}. 
Specially Fxie) is the orbital counting function of F. Let cp^ £ C^{G) be a 
nonnegative function supported on a small neighborhood of the identity 
element with /^^e = 1. Define a function defined on F\G (i.e., invariant 
under left F-action) by 

$,(F<7) = 5^0.(75). 
One is interested to estimate 



Jg 



IG 

'il)k{nxay)dny^^^^ dydk 

'M\K J y>T-^\\vok\\ J n^&{TnN M)\N M /M 

where ipk is an average of $e over a compact group M d K va Langland 
decomposition and {Tr\NM)\NM/M denotes the projection into T\G/M. 
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Geometrically {Tf^NA^)\NM/M denotes the quotient of a horosphere based 
at oo by the action of F n NM. See Sections [2] and fTOl for details. 

So it is important to estimate the average of a function over a horosphere 
at height < y < 1, /„^e(rn7VM)\iVM/Af V'('^xaj/)dn, for G L'^{T\G/M). 
Here one uses the matrix coefficient technique to draw Theorem 19.21 



for il) G L^(r\G)^, from which one deduces that Fx{e) ~ . The large 
part of the paper is devoted to justify this. 

First using a spectral gap, one can decompose L^(r\G)^ into Vs where 
the bottom eigenfunction (pQ is a unique ET-invariant function, and the rest 
V where V does not contain any complementary series Vs, s > sr- So any 
(j) e C^{T\G)^ n L'^{T\G/M) can be written as 



and hence {ay(j)i,4>2) = {4>i,<i)Q){ay4>Q,4>2) + Oiy^'''^ S2m{4>i)S2m{4>2))- See 
Corollary [521 

Secondly, we fix a positive t] G C^((r n NM)\NM/M) with r? = 1 on 
a bounded open set of (F n NM)\NM/M and vanishes outside a small 
neighborhood of B. Also for each e < eo, let be a nonnegative smooth 
function on AN~ M whose support is contained in 



and /^^ r^dv = 1. Finally define p^^^ on T\G which vanishes outside supp{'q)Uf 
and for g = nxan~m G supp{rf)We, 



This prj^e is introduced as a cut-off function to estimate the average of any 
function </> in C^(F\G)^ over a horosphere at height y in terms of e and an 
inner product of Qycf) and pr^^^- See Propositions 16.31 and 16.41 Finally in 
Theorem 19.21 one iterates the process to obtain the dominant term as the 
average of (pQ over the horosphere at height y, which is y^~^ ■ 

As an application to an Apollonian sphere packing in R^, we have 

Corollary 1.2. Given a bounded Apollonian sphere packing V in M^, the 
number of spheres whose curvatures are less than T grows asymptotically 
~ cT° for some positive a which is the HausdorfJ dimension of the residual 
set. 

Specially for asymptotic growth of the number of A;-mutually tangent 
spheres in with prime curvatures, we obtain: 

Theorem 1.3. Given a bounded primitive integral Apollonian sphere pack- 
ing V in M^, if ir^{T) denotes the number of k -mutually tangent spheres 




n-s 



4> = {4>,4>o)4>o + 4>' 



= {u, n A){u,, n N-){u,, n m) 



Pn,e{g) = r]{nx)re{an m). 
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whose curvatures are prime numbers less than T, then 

ttHT) « 



(iogr)'= 

for k < 5. 

Some of the lower bound seems to be known by experts like Sarnak 
for Apollonian circle packing. One can attempt a lower bound for r-almost 
prime curvature spheres as follows. The following corollary is a counting 
of 5 spheres kissing each other with certain properties, not the counting of 
individual sphere with prime curvature. See the last section. 

Corollary 1.4. Given a bounded primitive integral Apollonian sphere pack- 
ing V in , let iT^{Ty' denote the number of 5 spheres kissing each other 
(i.e. the number of orbits) among whose at least k curvatures are r-almost 
primes and all of whose curvatures are less than T where r is a fixed positive 
integer depending only on Apollonian packing. For any k < 5, 

For abundant literatures of this subject, see \16\ |17j for example. 

Organization of the paper: In sections 2 and 3, we give preliminary 
backgrounds for rank one symmetric space and specially introduce horo- 
spherical coordinates to calculate the Buseman function explicitly. In sec- 
tion [6l we deal with the bottom eigenfunction and its average over horo- 
spheres. This part is essential for our orbital counting problem. In section 
[71 we recall unitary representations of rank one semisimple Lie groups and 
generalize Shalom's result on matrix coefficents of spherical unitary repre- 
sentations of rank one group. In section [8] we give some consequences of 
unitary representations using spectral gap theorem due to Hamenstadt. Af- 
ter section [U we closely follow the proofs of [23] . We record them for the 
reader's convenience. In the final section, we are concerned with an inte- 
gral Apollonian sphere packing, and derive some asymptotic growth on the 
number of prime curvature spheres using uniform spectral gap theorem due 
to Bourgain-Gamburd-Sarnak, Varju-Salehi Gosefidy, Breuillard-Green-Tao 
and Pyber-Szabo. 



2. Preliminaries 

The symmetric spaces of R-rank one of non-compact type are the hy- 
perbolic spaces Hf, where F is either the real numbers M, or the complex 
numbers C, or the quaternions H, or the Cayley numbers O; in the last 
case n = 2. They are respectively called as real, complex, quaternionic and 
octonionic hyperbolic spaces (the latter one Hq is also known as the Cayley 
hyperbolic plane). Algebraically these spaces can be described as the corre- 
sponding quotients: SO{n,l)/ SO{n), SU{n,l)/ SU{n), Sp{n,l)/Sp{n) and 
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/ Spin (9) where the latter group of automorphisms of the Cayley 
plane is the real form of F4 of rank one. 

Following Mostow [28] and using the standard involution (conjugation) in 
¥, z ^ z, one can define projective models of the hyperbolic spaces as 
the set of negative lines in the Hermitian vector space F"''^, with Hermitian 
structure given by the indefinite (n, l)-form 

{z, w) = ZiWi H h ZnWn - Zn+lWn+1 ■ 

However, it does not work for the Cayley plane since O is non-associative, 
and one should use a Jordan algebra of 3 x 3 Hermitian matrices with entries 
from O whose group of automorphisms is -F4, see [28j. 

Another models of use the so called horospherical coordinates based 
on foliations of by horospheres centered at a fixed point 00 at infinity 
dH^ which is homeomorphic to (ndim^F — l)-dimensional sphere. Such 
a horosphere can be identified with the nilpotent group in the Iwasawa 
decomposition KAN of the automorphism group of H^. The nilpotent 
group N can be identified with the product F"^^ x ImF (see [28]) equipped 
with the operations: 

{tv)-{e,v') = {C + e,v + v'-2lm{{tO)) and {^,vy^ = {-^, -y) , 

where ((,)) is the standard Hermitian product in F""-*^, {{z,w)) = ^ Zjl/IJ. 
The group is a 2-step nilpotent Carnot group with center {0} x ImF C 
ipn-i jj^jp^ ^j^j g^pi^g itself by the right translations Th{g) = g ■ h , h, g S 
N. 

Now we may identify 

H^UdH^\{oo} -^N X [0,00) =F"-^ X ImF X [0,oo), 

and call this identification the "upper half-space model" for with the nat- 
ural horospherical coordinates (^, v, u). In these coordinates, the above right 
action of A^ on itself extends to an isometric action (Carnot translations) on 
the F-hyperbolic space in the following form: 

%),»^o)-(^'^'^) ' — > i^o + ^,vo + v - 2lm{{^,^o)) ,u) , 

where {^,v,u) G F"~^ x ImF x [0,oo). Hence rank one space can be 
written in horospherical coordinates as 

= {{z,t,y)\z e F"-\t E ImF,0 < y G M}. 

An hyperbolic isometry fixing and 00 acts as 

{z,t,y) {rzOw,r'^t,r^y) 

where Of is an element of M. Here the group MA fixing the origin of 
A^ and 00 is {U{n - 1) • [/(I)) x M in H^, {Sp{n - 1) • Sp{l)) x M in H^, 
and Spin{7) x M in Hq. Note that once the height y is fixed, which is a 
horosphere of H^, NM leaves invariant the horosphere, i.e., NM acts on 
each orbit of A^ by 

{O^,ix,v,0)) : iz,t,y) ^ {{x + z)Ow,v + t - 21m{{t,x)) ,y)- 
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Later when we use the notation (F n NM)\NAI/M, it denotes the pro- 
jection into T\G/M. Since G/M is the unit tangent bundle T^{G/K), 
(r n NM)\NM/M is the quotient of a horosphere under the action of 

r n NM. 

It is not difficult to show that {dt — 21m{{z, dz)))"^ + {{dz, dz)) is invariant 
under NM where ( ( , ) ) is the standard positive definite Hermitian product 
and ANM is a Borel subgroup fixing oo. Since a metric on can be 
written as Qy © dy^ where Qy is defined on a horosphere along a geodesic 
ending oo, one can give a metric [21] 

dy'^ + [dt - 2lYa{{z,dz))f +4y{{dz,dz)) 



ds' 



y2 



This metric has the sectional curvature between —1 and — | so that the vol- 
ume entropy for real hyperbolic n-manifold is (n — l)/2, the volume entropy 
of complex hyperbolic n-manifold is n and the volume entropy of quater- 
nionic hyperbolic manifold is 2n + 1 . 

The volume form on can be written as 
2n-i 

(1) y{n+i)/2 dVolzdy (real hyperbohc) 

(2) — —rdVolzdtdy (complex hyperbolic), 

yTl-\- V 

16""^ 

(3) 2n+2 dVolzdtdy (quaternionic hyperbolic) 

where dVolz is a volume form on F"^"*^. The set {x = {z,t)} is identified 
with the Nilpotent group in Iwasawa decomposition ANK of Iso{H^) 
and its right action is 

{z, t, y){w, s) = {z + w,t + s - 2Im((z, w)),y). 

The ideal boundary at infinity of is U oo. As usual A will be an one- 
dimensional group translating along {(0,0, y)|y > 0} and its right action is 
given by (z, t, y)ar = {rz, r'^t, r'^y) so that any point (z, t, y) G is 

(0,0,l)(^z,it)a^. 

In this way we will identify a point {z,t,y) in with [{-^z, ^t),a ^] G 
A^ X yl in a fixed Iwasawa decomposition KNA. 

Note that the metric gy © dy"^ is such that gy = e^^g^ © e^'^g^ on F"~^ 
part and ImF part, respectively, e~'^'^gy converges to nonriemannian metric, 
Carnot-Caratheodory metric on A^. It's distance is given by 

dN{{x,t),{w,s)) = \{x,t){w,s)-^\ = (|x-«;|^ + (t-s + 2Im((x,u;)))2)i/^ 

One can even define A^-invariant metric on horospherical model by 

d((x, t, 2/), {w, s, z)) = i\x -w\'^ + {t-s + 21m{{x, w))f + \y - z\^)^/^. 



EQUIDISTRIBUTION 



7 



3. BUSEMANN FUNCTION AND HOROSPHERE 

In this section we normalize the metric so that the sectional curvature 
is between —1 and — | and we fix a reference point o = (0,0,1). can 
be realized a unit ball in F". Two points (0',— 1) and (0', 1) will play 
a special role. There is a natural map from B"- to P(F"''^), where F"'^ is 
equipped with 

{z, w) = ZiWi H h ZnWn - Zn+lW^+l ■ 

defined as 

{W ,Wn) {w',Wn,l)- 

Prom B" to the horospherical model H", one define the coordinates change 
as 

/ / X . r ^' 2Im2:„ 1 - \zn\^ - \z'\^ ^ 
It's inverse from the horospherical model to PF"'^ is given by 

where v is pure imaginary, i.e., iv in complex case, and ivi + jv2 + kvs in 
quaternionic case. According to this coordinate change, (0', 1) = [(0', 1, 1)] 
corresponds to the identity element (0',0) in Heisenberg group, (0',— 1) = 
[(0',-l,l)] to oo, and (0',0) to (0,0,1). 

In the rest of the section, we carry out the calculations only in complex 
hyperbolic space but it goes through the quaternionic case. li x,y G 
and X,Y E C"'^ correspond to x,y, the distance between them is 

cosh^f^) - ^^'^^^^'^^ 
^ 2 {X,X){Y,Y)' 

where {X,Y) = Y17=i ~ Xn+iVn+i- A Busemann function based at is 
defined as 

B^{z)= lim {d{z,jt)-t) 

where 74 is a unit speed geodesic starting from o and ending ^ e dH^. First 
we calculate the Busemann function based at 00 = (0', — 1, 1). 7^ can be 
chosen as 

^(t) = (0', - tanh(t/2), 1) C C"'^ 
Then a straightforward calculation shows that 

d(.,7.)^i + log3%^, 

where z = {z',Zn) € B" and {{z,w)) = ^ZiWi is the standard positive 
definite Hcrmitian product. Wc denote it as a double bracket whereas ( , ) 
is (n, 1) Hermitian product on C"'^. Then 

^-("^=^"g l-(|.,'i)) - 



If we denote Q 
(4) 



(0', 
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-1,1) and Z - ^~ 1^ - '^"•1 
^-Bq{z) _ _ 



z, 1) in C"'^, it is easy to show that 



(Z,Q)(Q,Z)- 

Using elliptic isometry fixing o = (0,0,1) G "H", one can send Q to any 
other point in the ideal boundary. Since it preserves Busemann function 
and Hermitian inner product, above formula holds for any Z, Q. Using this 
formula, it is easy to check that e~^°°^^^ = y for z = {^,v,y) in horospher- 
ical coordinates. This explains the last height coordinate in horospherical 
coordinates. 

We are interested in e~^'^^^^ for Q ^ oo. Let Q = {(,,v) and z = 
{x,t,y) S Ti'^. These will correspond in C"'^ to {x, iH£j_z^±!i^ i+\x\_+y-tt^ 

— ) respectively. A direct calculation using equation 
4y 



and (e,i^,i±^ 
(HI) gives 



|2((x,0) 
Using the relation {{x — S,,x — ^)) 
formula becomes 

(5) e-^«(") 



X 



\x-e 

Ay 



— y + i{t — t>)P 

|x|2 + |^|2-2Re((2;,0), above 



(|x - + yf + {t-v + 2Im((x, ' 
Note that the denominator is comparable with the distance defined in section 

m 

4. Opposite Nilpotent group 

We fixed a Iwasawa decomposition ANK so that is a 2-step nilpotent 
group which is a Heisenberg group. To describe an opposite Nilpotent group 
N~ for later use, we introduce another coordinates changes. To define a 
Hermitian form of signature (n, 1) one can equally use 

1 

/„ 

1 

for the product. Somehow this matrix simplifies the calculations the most. 
In this context, a coordinate change (p from the Horospherical model T-L^ to 
p(F"'i) is given by 

(P{C,v,u) = [{-\C\'' -u + v)/2X,l\ 



see [2T]. 

One can easily show that the right translation by (t, t), i.e., the map 
(z, s) —7- {z,s){T,t) corresponds to matrix multiplication 



Z 1 



1 0' 

-T* / 

-H^ r 1 
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Since we consider right actions instead of left actions, we can define the 
opposite Nilpotent group N~ as 



N- = { 



1 0" 

T* I 
T 1 



-\r\^+t 



} 



2 

whereas consists of upper triangular matrices. Then using the multipli- 
cation rules in section [2] one can show that 

"^x^y — ^y'^x/yi^x^y — ^y^yx' 

This will be used in section [6l In these coordinates, the origin of the Heisen- 
berg group corresponds to (0,0, 1) in P(F"''^) and it is the stabilized by the 
right action of NM. This fact will be used in section [101 

5. Digression to a general Riemannian geometry 
Let M be a Riemannian manifold, u a function such that |Vn| = 1. Then 
Proposition 5.1. 

Au{x) = mean curvature at x of u~^{u{x)). 

If u is a Bus emann function on the rank one symmetric space X, u~^{u{x)) 
is a horosphere and Au{x) = D where D is the Hausdorff dimension of dX. 

Proof. 

An = divVu = div^, ^ = Vu. 

Since |C| = 1, {Vv^,0 = for any v G TM. Also since G End{TM) is 
symmetric (the antisymmetric part of is d{du)), 

hence (V^^,w) = = for any w G TM. So V^^ = 0. 

Let Ci be orthonormal basis of T^M with ei = ^. Then 



n 

divi = - ^(Ve^C, Ci) = - ^(Ve,C, Ci) = - ^ //(ei, Ci) 
1=1 i=2 1=2 

where // is a second fundamental form on u~^{u{x)). The second claim is 
well-known [T] (page 638-639). □ 

Note that for / : R M and n : M M, using div{gV) = {Vg, V)+gdivV 
for a smooth function g and a vector field V on M, 

(6) A(/(n)) = /"(n)|Vn|2 + /'(n)An, 

and if |Vn| = 1, A(/(u)) = f"{u) + f'{u)Au. 

For example if X is a symmetric space of rank one and n is a distance to 
a fixed hyperplane. An = (j){u) where is a solution to 

(piu) = Tr{II{u)) 
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where // : M — )• (n — 1) x (n — 1) symmetric matrices satisfying II' + 11^ + 
Riem = 0. Here TX = ¥Vu © (FVn)-^ and {Vu)^ = (/mF)Vn © (FVn)-^. 
Riem preserves this decomposition and 

Riem = —Aid on {ImF)Vu, 

= -Id on (FVn)-^. 
The solution of II' + 11^ + Riem = is 

(l)l{u)Id(^jjn¥)\'u + 4'2iu)Id(f'Vu)^ 

where + _ 4 = and + - 1 = 0. Then (^i(u) = ^^J^^u) ' '^2(w) = 
tanh(u). 

In conclusion 

2 

An = (/)(n) = Tr{II{u)) = {dim¥ - 1) — — - + {dim¥){n - 1) tanh(n). 

tanh(2n) 

For real hyperbolic 3-dimensional case, if ti is a distance to i/^, Au = 
2tanh(n). So if is an average of the bottom eigenfunction over H"^ at the 
distance u, 

A(p = (l}"{u) + (/>'(n)2tanh(u) 

and since 

= -5(2 - 5)<^ 

0(n) grows asymptotically 

ce-(2-5)« + de-^". 

If It is a Busemann function on the symmetric space X of rank one, then 
by equation ([6]) and An = D 

Aifiu)) = f"iu) + Df'iu), 

where D is a Hausdorff dimension of dX. Note that the solution space of 
the above differential equation is generated by e~™ and e,-^^-^)^ for some c. 
Hence if (jy^ is an average of the bottom eigenfunction (pQ on a geometrically 
finite manifold V\X over a horosphere at the distance n, 

<^^(n) = c<^„e-(^-^)"+d<^„e-'5^ 

But in later sections, we will use —A as the Laplace operator so that 

<(n) = c^„e(^-^)" + d^,e^". 

In section O we show that C(^q / 0. If Fq C F is a subgroup of index 
n, then l/-^/n0o is a unit norm eigenfunction of —A with eigenvalue 
5{D — (5) where 4)q denote the lift to Fo\X of on F\X. Then the average 
of 1/ y/n(j)Q over the horosphere at the distance n is 

%o^(n). 
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Specially the constant c^fq is equal to -^===0^^. This will be used in 
Section [TTl 

6. Bottom eigenfunction c/iq 

Let X be a rank one symmetric space and fix a origin (0, 0, 1) = o ^ X. 
Let T\X be geometrically finite with a critical exponent 8 > D/2 where D 
is the Hausdorff dimension of dX. Let be a Busemann function based at 
C, normalized that B(^{o) = 0. The Patterson-Sullivan measures 9x, x £ X 
satisfies: 

The measures Ox and 9y are mutually absolutely continuous and 
and for any 7 G F 

(8) 7*^x. = e^x. 

The function defined by 

Mx,t,y)= [ e-'^^^-^^'y^e 

J At 

where d9 = dOo is a fixed Patterson-Sullivan measure associated to a fixed 
reference point o S X, descends to a positive L^-function on T\X whose 
eigenvalue with respect to the Laplace operator — A is 5{D — 5) where D is 
the Hausdorff dimension of dX. We always normalize it that its L'^-norm 
is 1. The function (/iq is given by in horospherical coordinates according to 
equation ([5]) 

<l)oix,t,y) = [ {- — .2 I 1 oT // A\\\2 )'^^^(^'^) 

Ar (f - Cr + yy + it-v + 2Im((x, C)))^ 

Note that the formula for real hyperbolic case seems a bit different from this 
one but it is due to the fact that the C sits as a Klein model [15| (not 
Poincare model) and that the curvature is —1/4 in this section. After we 
normalize the metric back to between —4 and —1, all the formulas will turn 
outright. See section [71 
Let 

fPo iv) = i (poix,t,y)dn 

J{rnNM)\NM/M 

be the average over a horophere where dn = dtdx. This is independent of 
the choice of a fundamental domain since dn is A^M-invariant and 0o is 
P- invariant. 

We will be interested in horospheres whose images are closed in T\H^. 
Specially we will consider horospheres whose base points are either parabolic 
fixed points or points outside the limit set of P, see Lemma fl 0.11 
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Proposition 6.1. We identify N with a horosphere based at oo (in fact, 
the orbit of N) and suppose (T D NM)\NM/M, which is a quotient of 
the horosphere by F n NM , has a closed image in T\X. Then (j)^ (y) ^ 
y^~^ for all Q < y <^ 1 where D is the Hausdorff dimension of dX, i.e., 
(n — l)/2,n,2n + 1 for real, complex and quaternionic hyperbolic space re- 
spectively. 

Proof. We carry out the calculation in complex hyperbolic case and indicate 

the difference in the other two cases. 

CASE I) If oo ^ Ar, then NM n T is trivial and 

4>o{y)= I / ( I ,12^ \ Yde{c,v)dn. 

r V y ; -T V y ) 

Change the variables to 

^Jy y 

to get 

dtdx = y^dsdz 

and 



-^,s = -, X e C""\t G 



(9) <P^{y) = A'y^~' 

f dsdz 

^ ^ Ar [(k - c? + If + \s-v' + 2im{{z,cm' ■ 

Change {w,t) = {z, s){C' ,v')~^ . Then since dn is A^-invariant 



dtdw 



In [\w\^ + \t\^ + 2\w\^ + 1Y ' 

The equality holds since d9 is a probability measure. One can show that 
the integral converges for 5 > ^ to a nonzero number. In this calculation, 
note that for real hyperbolic case, there is no t factor, so from z = x/^, 

dx = {,yy)'^~^dz, this gives y~ ^ in front. In quaternionic case, dx = 
(y/y)"^""^, dt = y^ds, to get y2n+i-<5 -^^ front. CASE II below is similar. 
CASE II) oo G Ar 

By the theorem of [B], oo is a bounded parabolic fixed point, so r\(Ar — oo) 
is a compact set. Hence it is bounded. Let F\ = {T r\ NM)\Ar C F/v = 
(r n NM)\NM/M be fundamental sets under the left action of T n NM. 
Note that for a fixed Fa C and X G , ( G Fa, by the property ([8]) of 
Patterson- Sullivan measure, 

{n:Hxm = en-^x{0 = ex{nQ) 
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for any n G NMr\T. Also by the property ([7]) of Patterson-Suhivan measure 



dOo 

dO, 



dOn 



This imphes that for Q G i*A 



uFa JFn J Fa Jti-^Fn 



Hence, 



J Fai J Ay nernNM "-^A ■J Fn 

^^T.rJFA Ju-^Fn JfaJn 



In terms of the coordinates 

(11) ^^{y) = 

[ ^ 

lirnNM)\Ar [{\z-C'\^ + ir + \s-v' + 2lm{{zX'm' ' 



N 



where {C,v') G aj_F\. A similar estimates holds to conclude that 

Vv 



□ 



Fix generators fi, • • • ,ffc in NM/M corresponding to the axes of screw 
motions of A^M n F so that vi,--- ,Vk together with Vk+i,-" ,V2n-i are 
basis of NM/M. Denote A^-*" the subspace generated by Vk+i,--- ,V2n-i- 
Let F\ d B C F/v be an open set such that 

(1) if oo ^ Ar then eo(-B) = mlu,^FA,x&B'= \x - u\n > 

(2) if oo a bounded parabolic fixed point then eo(-B) = miu£FA,x£B'= \x — 
u\]^± > where B" = \ B. 

For such an open set B 

Proposition 6.2. If6>§, 4>^ (y) = Mx,t,y)dn + O^^^B)iy^) and 

4>^iy) = 0{y^-'). 

Proof. We give a proof in complex hyperbolic case but as in the previous 
Proposition, the other cases are similar. When cxd ^ Ap, using equation Q, 
we want to estimate 

4<5 n-5 f f d,sdz 

^ JBJAr[i\^-c? + ir + \s-v'+2im{{z,c'm' ■ 

Since 

1 



[i\z - CP + 1)2 + \s-v' + 2Im((z,C'))l 



/\\|21<5 
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- [\z-CT + \s-v' + 21m{{z,C')W]' [dN{iz,s),{C',v'))^]^' 
Since dn is invariant under A^, the above integral is 



where w = {z, s){C' ,v')~^ in Heisenberg group. Let w = {x,t) and write 
dtdx = r'^^~'^drdtdS where dS is a volume form on unit sphere in M^""^. 
Then 

^Syn-S [ dn ^ ^Syu-s f r^--HrdtdS 



r'^^~^drdt 



<Cy 



n—5 



r>tn 



5 ■ 



Letting t = tan^r^, and for 6 > ^ it becomes 

n-5 f r^^'-'dr f de 



As before we have 

.5..n-5 r r^'^-^drdtdS 



Cy' 



r,t>o {r' + t^ + ly 
r'^'^^^drdt 



r>Jt>o (r4 + 2r2+t2 + i)5- 

Hence if 4(5-2-2n+3 > 1, (i.e. 5 > f ) then by letting t = tan 6'v^(r4 + 2r2 + 1) 
the integrals converge to a nonzero number to conclude that 

When cxD is a bounded parabolic fixed point, the similar estimates holds. 
Using 



= ^ / / e-^^o^'^'^'yUndO = / / e-^^o^'^'^'yUnde, 
and following the equation pT]) . we want to estimate 

aSu-S f f ^ ,n 

^ L,rn.MnBJ(rnNM)\Ar [{\z - + + \s - V + 21m{{z , (WY ' 

where {C',v') E aj_F\. The same estimation as in oo ^ Ap gives 

Vv 

f e-^^o^^'^'^yUedn < 

B'' JAr 
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to get 



{y)= f [ e-^^o'''''*'yU9dn+ [ [ e-^^o'^'^'^'yUOdn 
JB JAr JB" JAr 

J B JAr 

Also in this case the similar estimates give that for S > ^ 

□ 

We fixed Iwasawa decomposition ANK so that is a Heisenberg group. 
Let be the opposite Nilpotent group to N so that 

NxAyiN-yiM^G 

is a diffeomorphism around a neighborhood of e and dv is a smooth measure 
on AN~ M so that dn ® dv \s a. Haar measure d^i on G. Fix a left invariant 
metric do on G and is an e-neighborhood of e in G. Since AxN xK ^ G 
is a diffeomorphism and hence a bi-Lipschitz map around the neighborhood 
of e, there exists / > such that is contained in Ai^Ni^Ki^ once we fix 
some eo and take e < eo- We fix a positive r/ G G'^{{N M r\V)\N M /M) with 
r/ = 1 on a bounded open set B of F/v so that eo(-B) > as in the previous 
Proposition 16.21 and vanishes outside a small neighborhood of B so that 

<^o{y)= I Mri,yHn)dn + 0{y'). 

JB 

Shirinking eo if necessary, we further assume that 

supp{r]) X {Ueo n AN~M) T\G 

is a bijection to its image. 

For each e < eo, let be a nonnegative smooth function on AN~ M whose 
support is contained in 

= {u, n A){u,, n N-){u,, n m) 

and Jyy r^dv = 1. Finally define on r\G which vanishes outside supp{r])U^Q 
and for g = nxan~m £ supp{rj)Wt^, 

PvAa) = ri{nx)r^{an~m). 

Then 

Proposition 6.3. For small e ^ eo and for y < 1, regarding (po G C(r\G)^ 

cP^iy) = {aycPo,Pr,,e)L'^+Oiey''-') + 0{y'). 
Proof. As usual we give a proof only in complex hyperbolic case. Since 



{ay4>Q, Pr),e) L'^ = I r^{h)l (j)o{nhay)ri{n)dndu{h) 

JW^ J {NMr\T)\NM/M 
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= / (f>o{nhay)ri{n)dn 
J{NMnr)\NM/]\i 

we need to estimate (j)o{nhay) for n G Fn and h G We- 

For h = ay^n~m S so that |yo ~ 1| = 0{e), 

Since n~^. = ay^Ux^ki G Aiy^Niy^K so that |yi — 1| = 0{ye), 

flhciy — ^^xiyyoyi^yyoyi^l'^' 

Since (po is iT-invariant and dn is A^-invariant, 

(l)o{nhay)rj{n)dn = / (t)Q{nnx^yyoy^ayy^^y^)ri{n)dn 
Fn JFn 

4'o{nayyoy^){r]{n) + Or^{e))dn 

Fn 

by letting n' = nn^iyyoyi and so ??(n) = ri{n' {n^^yy^^y^Y'^) and as |r/(n') 
rii'lT'' (nxiyyoyi = 0^(e). By Proposit lon lD.zl 

(t)o{nhay)ri{n)dn = / (j)o{nayy^-^y^)r|{n)dn + 0^(6(1)^ {ayy^y^)) 
Fn JFn 

= (t>o (yyoyi) + 0{{yyoyiY) + Or^{e(j)Q {ayyoy^)). 
Since by Proposition 16.21 

^.^ t/*^ (yypyi) ^ ^.^ '/'o' (yyoyi) {yyoyi)^~^ y^~^ 



y-^'o ^oiy) y^o{yyoy^)D-S yD-5 

D-5 

nhay)rj{n)dn = + 0(/) + O^iey""-^). 

Fn 



limiyoyif-' = {l + 0{e)), 

y-^O 

and since c/)^ [ayyoy^) = O((yyoyi)^"^) we get 



□ 



Let {Zi, • • • ,Zk} denote an orthonormal basis of the Lie algebra of G and 
r C G a discrete subgroup. For / G C°°{T\G)^ n L'^iV\G), one considers 
the Sobolev norm Sm{f)'- 

SmU) = max{||Z,, • • • Zi^{f)\\ : 1 < i,- < k]. 

The following is standard. For (j) G C^{T\G)^ , there exists G C^{T\G)^ 
so that 

(1) for small e > and h G [7^, 

10(5) - -^(5/^)1 <e0' (5) 

for any g G r\G. 

(2) by Sobolev embedding theorem, there exists q so that Sm{(p') ^ 
Sq{(j)) for each m, where the implied constant depends only on supp(0). 
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Proposition 6.4. Let cj) G C°°(r\G)^. Then for any < y < 1 and any 
small e > 0, 

\Iv{4>){ay) - I < {e + y)Iri{(t)'){ay) 

where /^((/))(ay) = / (j){nay)ri{n)dn and i] £ Cc{{NM n T)\NM/M). 

Proof. For /i = an~m G W^:, an~may = aya{ay-\n~ ay)m. Hence nan~may = 
nay{aay-in~ Oym) where Oy-in^ay € Uy^^ D . Then as (j) is if -invariant, 

\4>{nay) - 4>{nhay)\ = \(t){nay) - (t){nayh')\ < <p'{nay){e + yeo) 

where h' = aOy-m^ay G {U^ n A){UyeQ n N^). Hence 

\(j){nay) - / (t){nhay)re{h)dv{h)\ <.(t){nay){e + yeo). 

By integrating on {NM n T)\NM/M we obtain 

\Ir^{4'){ay) - {ay(j3,pr^^^)L2(r\G)\ < (e + yeo)/,,((/>')K)- 
Since eo is fixed, we get the desired result. □ 

7. Unitary representations of rank one semisimple group 

From this section, we normahze the metric so that its sectional curvature 
is between —4 and —1. Equivalently we have to multiply the metric ten- 
sors ds^ by 1/4. Then the volume form for complex hyperbolic n-space is 
multiplied by 2~^"'. Then the distance from (0,0, 1) to (0,0, y) will become 
log y^. So by changing the variable = w and abusing the notation by 
putting w back to y, the volume form for complex hyperbolic n-space is 

(13) 

For real hyperbolic n-space, 

(14) \dVohdy 

yn 

For quaternionic hyperbolic n-space 

1 

8y^ 

where dVolz is a volume form on F""^. Also in all the formulas in section 
[U n should be read as 2n under this normalization. 
To O] 

Let T] = 1 be in 5o(n, l),5u(n, l),sp(n, 1) so that e^^ is the 1- 
|_0 1 oj 

parameter subgroup constituting A in KAN. A direct calculation shows 
that 

Qk = Ker{adr) - kid), k = 0, ±1, ±2 

are only root spaces and there are only two positive roots /3, 2/? (in real 
hyperbolic case, 2/3 is not a root) so that /?(??) = 1. Note that real dimension 



(15) -^dVohdtdy 
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of Qi is n — 1, 2n — 2, 4n — 4 resp. and that of 02 = 0,1,3 resp. Then the 
half sum p of the positive roots is 

Since we wih work with horospheres based at oo and expanding ones as 
y — 7- 0, we win define a positive Weyl chamber by 

A+ = {ay\0 <y<l} 

as a multiphcative group. So — /?, —2/3 will be positive roots in this paper, 
which will change the plus sign to the minus in all formulas in the literature. 

In this section, we prove the following. The kind of estimate we look 
for was first established by Cowling, Haagerup and Howe for tempered rep- 
resentations (i.e. unitary representations weakly contained in the regular 
representation). A unitary representation (F, vr = J^^^TTxdn{x)) of G is 
called tempered if one of the following is true, see for example [9| l29]. 

(1) For any /C- finite unit vectors v and w{i.e., the dimension of the 
subspaces spanned by Kv and Kw is finite), 

\{7T{g)v,w)\ < {dim{Kv)dim{Kw)y/^EG{9) 

for any g £ G where He is the Harish-Chandra function on G. 

(2) For almost all x € G, the irreducible representation tt^ is strongly 
L?''^'^ {G / Z {G)) , i.e., for any p > 2, there exists a dense subset W C V 
such that for any v,w ^ W, the matrix coefficient g — )• {'7Tx{g)v,w) 
lies in LP{G/Z{G)). 

(3) For almost all x £ G, tTx is tempered in the sense of (1). 

(4) TT is weakly contained in the regular representation L^(G), i.e., any 
diagonal matrix coefficients of vr can be approximated, uniformly on 
compact sets, by convex combinations of diagonal matrix coefficients 
of the regular representation L^(G). 

The following theorem is due to Y. Shalom, [35], Theorem 2.1 p. 125 in 
case G = SO{n, 1) or SU{n, 1). We generalize it to other rank one groups. 
For the notational simplicity, we fix the notations first. 

Let G be a M-rank one simple Lie group. Pick an Iwasawa decomposition 
G = KAN. Let A G be a complex linear form on the Lie algebra of A. 
This gives rise to the character 

on A. Let M be the centralizer of A in K. Let Z'^ denote the space of 
i^-finite complex valued functions on G such that f{gman) = a~'^ f{g) for 
all a G A, m G M and n £ N (notation taken from [24]). G acts on Z'^ by 
{TTx{g)u){h) = u{g~^h), for g, h e G. 

Notation 1. Let p denote the half- sum of positive roots. We write A = 

(1 + S)p with S & C, Zg = Z^, TTg = 7T\. 
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A description of the set of s G C such that Zg is irreducible and admits 
a G-invariant inner product can be found in [24], Theorem 10 page 641. 
It is the union of the hne {Re(s) = 0}, and of the closed critical interval 
CI = [— si(G), a symmetric interval on the real line. It is shown in 

[23], that for < A < p, tta and Ti2p-\ are equivalent. The representations 
obtained for {Re(s) = 0} are called spherical principal series representa- 
tions of G. The representations obtained when s G CI are called spherical 
complementary series representations of G. 

The word spherical refers to the fact that these representations contain a 
unique iiT-invariant line, generated by the function vx which equals 1 on K. 

According to Harish-Chandra, [19], every irreducible unitary representa- 
tion of G which contains a nonzero i('-invariant vector is isomorphic to a 
representation from the spherical principal or complementary series (except 
possibly for the trivial representation). 

Let G be a M-rank one simple Lie group. Let (tTs, Z^), s G [0, si{G)], be a 
representation belonging to the spherical complementary series of G. Pick 
a unit A'- invariant vector Vs in Zg. Define the function 

^s{g) = {Trs{g)Vs,Vs)s- 

When s = 1, A = 2p, hence 7r2p is equivalent to vro. So Hi = Harish- 
Chandra function on G. 

If KA^K is a polar decomposition of SU{n, 1), then the Haar measure 
on it is ([35j) 

where log : G — t- g is the inverse map of the exponential map and a{g) is 
the component of in the polar decomposition of g. Note that in this 
formula, A'^ is regarded as an additive group. Then for y < 1, the measure 
is comparable to e~'^^^^^°^°'^ dkdadk, that if we use logy = a G A'^ to make 
A"*" a multiplicative group so that do = ^ it is 

(16) 

2n+i dkdydk. 

We hope that this switch between multiplicative and additive group does 
not cause any confusion to the reader. This is consistent with the volume 
form (jl3p on H'^. For real hyperbolic space, Haar measure is 

^g-/3(loga) _g/3(loga)^n-l^^^^^^ 

and comparable to 

for y < 1, hence after put logy = a, the Haar measure on KA^K is 

— dkdydk^ 
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which is comparable. In any case one can write the Haar measure on KA^K 
as 

(17) —^y^dkdydk 

for y < 1. For quaternionic hyperbohc case is the same. 

For A G ia*, and for log 0(5) < 0, the principal series ttx has matrix 
coefficient decaying rates for X-invariant unit vector v\, 

\{7rx{g)vx,vx)\ < Eig) < C\l - /3(loga(5)|ef /^C^^'^^^)). 

See dU section 4.6.4. For < A < §/3 G 0*, and for log a(c/) < 0, the 
complementary series representation has matrix coefficient decaying rates 

(18) < C7|l-/3(loga((7)|e(-^+f/3)(^°sa{9)). 
See [55] (equation (10) page 132) also. 

Proposition 7.1. For all g £ G, there exist c{G),C{G) such that 

(19) c{G)E,{gy < Es{g) < C{G){1 - logE,{g))E,{gr . 

Proof. Let A = (1 + s)/3, < s < si{G). Then tt\ is equivalent to n2p-\ 
where 2p — X = p — sp. In view of equation (fTSj) . 

Es{g) < C|l-/3(loga(5)|e(-2p+A+p)(ioga{g)) = c\l - ^{loga{g)\e^'P'>^^°^''^^^\ 
But it is known [14j (Theorem 4.6.5) that 

Hi(a) = H(a) >e^('°sa)^ a G A+. 

Hence we get 

Es{g) < C{G){1 -logEi{g))E,{gy. 

For the lower bound, in [14J (4.7 (4.7.13)), [35] (equation (11) page 132), 
it is shown that 

^/g(p-A)ioga(9) <Ex{g), < A < p. 

Hence again using the fact vr;^ is equivalent to 7r2p_A) we get the desired lower 
bound. □ 

Theorem 7.2. There exists a constant C{G) such that for all K-finite vec- 
tors u, V £ Zg = Z^, 

{7rs{g)u,v)s < G{G) (dim Span(ifn) dim Span(i^t;))^/^|n|s |t^|s Hs((7). 

Proof. The main ingredient in the proof is Cowling, Haagerup and Howe's 
temperedness criterion {9] . The following statement is a combination of their 
Theorems 1 and 2. 

Proposition 7.3. Let tt be a unitary representation of a semi-simple alge- 
braic group G. Then the following are equivalent. 

(1) vr has a dense set of vectors whose coefficients belong to L^+'^(G) for 
all e > 0. 
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(2) for all K-finite vectors u and v in tt, 
(20) {iT{g)u,v) < C(dimSpan(i^n)dimSpan(Kt;))i/2|^| \v\Ei{g) 

where Si = S is the Harish- Chandra function on G. 

Also, by Proposition 110. 2^ for all g £ G, 



Shalom's trick consists in tensoring representations until they become 
almost square integrable. Let s G [0,si(G)]. For G = Sp{n,l), n > 2, 
•^i(^) = liri > ^- "^^^^ ^^^^^ exists t G [0, si(G)] such that 1 - s - t G 



[0,si(G)]. For G = F^^^, si{G) = ^ > |. Thus there exists t e [0,si(G)] 



such that 1 — s — 2t E [0, si{G)]. Consider the unitary representation 

Vr = VTs TTj (g) 7ri_s_j (resp. tt = TT^ (g) TTt (g) TTf (g) 7ri_s_2t). 

Since all vTs are irreducible, the translates of Vs generate a dense subspace 
of Zg. Let til, vi belong to this subspace. Then {TTs{g)ui,vi)s = 0(Es{g)). 
Similarly, pick U2, V2 (resp. u^, v^, u^, V4) in the vector space generated 
by translates of vt (resp. of vi-s-t, vi-s-2t)- Let U = ui (g ii2 g) U3, V = 
f 1 <g W2 ^3 (resp. til g) U2 (g ^3 (g U4 etc.). Then 



{7T{g)U,V) = 0{Es{g)Et{g)E,^s-t{g)) (resp. 0{Es{g)Etig)Et{g)Ei_s-2t{g))) 
= 0((l-logHi(5))^Hi(5)). 



thus belongs to L^+''(G) for ah e > 0. 

Since these products generate a dense subspace of the tensor product. 
Proposition 17.31 applies, vr is tempered, and inequality (pOj) holds for all in- 
finite vectors in the tensor product. Let u, v £ Zg be i^-finite vectors. 
Since 

{'iT{g)u (g -ut (g vi-s-t,v (g ft (g vi_s-t) = {T^s{g)u, v)s'^t{g)'^i~s-t{g) 



(21) 



c{G)E^{gr < Eg{g) < C (G) {1 - logE,ig))E, {gY . 



Since y^logy — )• as y — )• for any 6 > 0, for ay,y < 1 we have (1 — 
log Hi((7))'^Hi((7) < GyP^^ for any 5 > 0. Since the Haar measure is given by 




and 



ti (g (g Vl-s-t\ = \U 



-U (g -Ut (g Vi_s_t\ = \v\s-, 



{TTs{g)u,v)s < C (dim Span(inn) dim Span(irti)) 



1/2 



Si (5) 



U\s \V\s — 



Et{g)Ei^s-t{g)' 




Ht(5)Hi_,_t(5) - c(G)2 - c(G)3' 
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yielding the announced inequality. □ 
So for Oj/, y < 1, from equation (fTSj) 

(22) \{ayVx,vx)\<C\l-\ogy\y^-^ 

We will change the parameter for A from 0<A</9toy/3 = p<A< Dfi = 
2p so that the bound becomes 

C|l-logy|y^-^ 

and when A = I? it represents a trivial representation and when X = D/2 
represents a principal series. Since logy — as y — for any e > 0, we 
will write 

(23) \{ayVx,vx)\<Cy''-^-^ 
for p < A < 2p. 

8. Bottom eigenspectrum of Laplace operator 

Henceforth, we use the notation (VxtTTx) to denote the spherical principal 
or complementary unitary representation of rank one group G. Let X be a 
rank one symmetric space and fix a origin (0, 0, 1) = o G X. Let M = T\X 
be geometrically finite with a critical exponent 6. Let Bg be a Busemann 
function based at 9 normalised that Bq[o) = 0. The function defined as 

U^,t,y)= I e-'^^^-''^y^de 

where d9 is a fixed Patterson-Sullivan measure associated to a fixed reference 
point o G X, descends to a positive function on M whose eigenvalue 
with respect to the Laplace operator is 6{D — 6) where D is the Hausdorff 
dimension of dX. For H^, H^, H^, D = n — l,2n,4n + 2 resp. under the 
normalization of the sectional curvature between —4 and —1. Note that 
L'^{r\G)^ is naturally isomorphic to L'^{r\H^) by averaging over if -orbits, 
and this isomorphism intertwines the action of C, the Casimir operator, 
with that of —A. For geometrically finite groups, the relation between the 
bottom spectrum Aq and the critical exponent 5(=Hausdorff dimension of 
the limit set) is 

Xo = 5{D-6). 

In this paper, we will assume that 6 > D/2 and {Vx,Tr\), D/2 < X < D 
to denote the spherical principal or complementary unitary representation 
of rank one group G. Fix D/2 < sr < S so that there is no eigenvalue of 
the Laplace operator between sr{D — sr) and the bottom spectrum 6{D — 
6) in L'^{T\G/K). Such sr exists since the spectrum is discrete for small 
eigenvalues, see [I9j for general geometrically finite manifolds with pinched 
negative curvature. Then 

L^{T\G)^ = V5®V 
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where is a complementary series corresponding to 5 and V does not 
contain any complementary series Vg for s > sr- 

If Vs is an irreducible factor for V, we want to estimate the matrix coef- 
ficient decaying rate. 

For Xi, an orthonormal basis of the Lie algebra of K with respect to an 
j4(i-invariant scalar product, let co = ^ — '^Xf. This is a differential operator 
in the center of the enveloping algebra of Lie(K) and acts as a scalar on each 
X-isotypic component of Vg- 

Proposition 8.1. Let (VjTt) be a unitary representation of G = SO{n, 1), 
SU{n,l), Sp{n,l) which do not weakly contain any complementary series 
representation Vg for s > sq. Then for any e > 0, there exists q such that 
for any smooth vectors wi,W2 G V, and y < 1, 

\{aywuw2)\ < c,y^-^«-^||a;™(u;i)||||L^'"(u;2)||, 
where Am > rank(K) + 2^{positive roots}. 

Proof, vr decomposes as Jq®^" Pzdv where G is the unitary dual of G, mz 
is the multiplicity of an irreducible representation pz, and is a spectral 
measure on G. Then pz is either tempered or isomorphic to a complementary 
series Vg for D/2 < s < sq < n - I, D/2 < s < sq < 2n, D/2 < s < sq < An. 
Note that by Kostant, sq cannot exceed An for quaternionic case, 
Then for if-finite unit vectors wi and W2 of vr, if we write 

Wi = ®v^, G Vx, 

then by Theorem 17.21 

{ayv^,v^) < C(G)| |t;^| I \\v^\\ (dim Span{Kv^) dim Span{Kv^))^/^Ex{ay). 

By equation ([23l) . using Cauchy-Schwarz inequality and dimSpanKwi > 
dimSpan(i^T;^), finally we get 

\{aywuw2)\ < C(G)y^-^°-^nVdim(Kw;i)- 
From i^- finite vector to smooth vectors, it is standard, see [25], that 

\{ayWi,W2)\ < C{G)y''~'»-^\\u"^{wi)\\\\u;^{w2)\\, 

where Am > rank(K) + 2^T,'^ and S+ is a set of positive roots of K. In 
real hyperbolic n space case {K = SO{n)), 

m = 0,1, 2, A 

for n = 2, 3, 4, 6 and for n > 5 (odd), m = {n — 1)^/4 and for n > 8 (even), 
m = /A. □ 

From this we get 

Corollary 8.2. Let V he a geometrically finite discrete subgroup of G with 5 
as in the standing assumption. Then for any 0i G G°°{T\G)^nL'^{T\G/M), 
(t)2 G C^{T\G/M) andO<y<l, 

{ay(j)i,(l)2) = {4>l,(po){ay(j)0,<l)2) + 0{y^~''^ S2m{(t>l)S2m{(t>2))- 
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Proof. Note first that 

l'^(t\g)^ = Vsev 

where is a complementary series corresponding to 6 and V does not 
contain any complementary series Vg for s > sr- Put (f>i = (t>o)4>o + 4>t- 
Since (/>o is K-invariant, is also K-invariant. Then 

{ay(l)i,(j)2) = ((/>!, <Ao>(«i/<Ao,<A2> + {ay(j)i,(l)2) 

since S2m{4>i~) ^ 'S'2m (</>«) and by Proposition 18.11 □ 

9. Equidistribution 

Suppose the image of a horosphere is closed in T\H^ a geometrically 
finite manifold. Let N{J) = {[n] £ {NMnT)\NM : T\TnAnJ / 0} where 
J C r\G a compact set, i.e., the set of elements on Fjy whose orbit under 
the A-flow intersects J. We use the notations of section [2] and section [6l 

Lemma 9.1. N{J) is bounded if J is compact. 

Proof. Let Cr be the convex hull of the limit set and r\Cr the convex core 
whose volume is finite due to the geometrical finiteness of F. We claim that 
there exist L ^ 1 and a fundamental domain i*" of F such that 

(24) {{z,t,v) : {z,t) G FN,\{z,t)\ > L ovv > L} C F. 

If the base point oo of the horosphere is not in the limit set of T, since the 
domain of discontinuity is open, there exists a horoball based at oo which 
embeds into the manifold T\Hf, hence embeds into a fundamental domain 

F of r. 

If the base point is a bounded parabolic fixed point and Fjy is unbounded, 
since F is geometrically finite, the intersection Cr H {Fjy x [0, cxd)) C {{z, t) G 
F/v : I (-2,^)1 < X [0, oo) for some K. Otherwise since a horoball H = 
f) : V > R} \s stabilized by F n NM for some large R and quotient 
down to a cusp of T\Hf, the quotient of H Cr would have an infinite 
volume, hence the volume of the convex core would be infinite. Take a F- 
invariant small e-neighborhood M of Cr. Then Mr\ {Fjy x [0, oo)) C {{z, t) G 
F/v : I (-2,^)1 < K'} X [0, oo) for some K' . The nearest point retraction 
r : {H^ \ M) — )• dM commutes with the action of F. Since dN is invariant 
under F and F/v x [0, oo)\AA C r~^(57\An(F/v x [0, oo)) due to the convexity of 
N , if we denote J- a fundamental domain of F on dN , and Fc a fundamental 
domain of F on Cr, 

r-\F)yjFc = F, 

and so for some large L S> 1 

{(^,t)GF^:|(z,t)|>L}x[0,oo)cF. 

When Fjsi is bounded, there is nothing to prove. 

Now suppose N{ J) is unbounded. Then there exist sequences rij G -F/v — 
oo,aj G A,^j G F and Wj G J such that rijOj = "fjWj. As J is bounded, we 
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may assume that wj ^ w £ G. Take 70 S F such that 7ou;(0, 0, 1) G intF by 
choosing the base point different from (0, 0, 1) if necessary. Then for large 
j, 70^^(0, 0, 1) G intF. Since rij — )• cxd, by (f24l) . njaj{0, 0,1) £ F for large j. 
Therefore for large j 

7o77^njaj(0, 0, 1) = -foWj{0, 0, 1) G intF n 7o77^(^)- 

Since F is a fundamental domain of F, 70 = 7j for all large j. Since Wj — >■ 
ly, njttj = 'jjWj is bounded, a contradiction. This shows that N{J) is 
bounded. □ 

So for (p £ Cc{T\G/M) and any function rj G Cc((A^MnF)\A^M/M) with 
v\N{supp<t>) = 1, we have 



4>{nay)ri{n)dn = / (f){nay)dn 
{NMnr)\NM/M Jn{supp4>) 

since 4>{nay) = for n outside N(supp(j)). Hence we may assume that 
I(NMnr)\NM/M (t>{'nay)dn = I,^{(f)){ay). 

Since is an e-approximation in ^-direction S2m{Pe) = Or^{e'''^) for some 
p > 0, where Pe = p^^^. Fix / an integer so that 

{D-6)ip + l) 
(S-sr) • 

Let V'o = and define 

V'i = V'^_i, l<i<l 

inductively. 

Then a similar technique in [23j shows that 

Theorem 9.2. For any (p G C^{T\G)^ , 

[ (t>{ny)dn = (0, 4>o)c^,y''-'{l + 0(/) 

for some constants C(f)f^ and 5' . Here 6' depends only on the spectral gap and 
the implied constant depends only on Sobolev norm Sm{<P) and rj, hence on 

Proof. Now for each < i < I — 1 

-^r,(V'i)K) = (o-yAiPe) + 0((e + y)I^{ipi+i){ay)) 
by Proposition 16.41 and 

^viMay)) = {ayi^i, p^) + 0{{e + y)Ir,{ij'i){ay)). 

Heve since Irj{ij'i){ay) = J^NMnT)\NM/M ^i('^^y)'li'^)'^'^ ^ <^ /(7VA/nr)\JVM/M ^0 
where C = max ip'i(nay) on supp{r]), we map write 

IviMo-y)) = {ay^uPe) + Oriiie + y)) 

where the implied constant in O^^ depends on f(^j^Mnr)\NM/M 'n{''^)dn. 
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By Corollary 18.21 for 1 < i < ^, using Sktpi <^ Sqipi^i, 

{ay^i,p^) = {llJi,(l)Q){ay(j)o,Pe) + Oiy^'"^ S2mitpi)S2m{Pe)) 

= 0{{aycPo,Pe)\m2) + 0{y^~^^e~PSg{cP)) 
= 0{{aycPo,Pe)S2m+im + 0{y''-'^e-PSg{cl>)). 
Hence for ?/ < e, 

/-I 

Irj{(l)){ay) = {ay(P,pe) + 0{{ay'il^k,pe){e + yf) + O^(/^(V'0(aj/)(e + v)^) 

k=l 

= {aycP, p,) + Oi{ay^o,Pe)eSgi<P)) + 0{ey''''^e-PS,m + 0,(e') 
= </.o)(aj,0o, Pe) + 0{{aycl)o,Pe)e) + 0(2/^-^^^-?) + Qie^) 
= (0, 0o)0o^(a,) + Oiy') + 0(6y^-^) + 0(y^-^r,-P) + Qie') 
by Proposition 16.31 All the implied constants depend on Sq{(j)) and / rjdn. 
Setting e = y('5-«r)/(p+i)^ have e' < y^'^ and 

</.(na,)(in = I,(</<)(a,) = (0, 0o)0^K)(l + 0(y(^-^i^)/(P+i))) 

(ArMnr)\AfAf/M 

using the fact that <pQ (ay) = c^^^y^"^ + df^^y^ and the fact 5 > D/2 > 
^ + - 5. The claim follows by setting 6' = □ 

There exists a Burger- Roblin measure p onT^X = G/M which descends 
to T^{T\X) and satisfies 

p[tP) = (V',0o) 

for ET-invariant functions Tp G Cc{G). Specially /i((/'o) = 1- Indeed one 
can write down jl explicitly as follows. A fixed Patterson-Sullivan measure 
on dX = K/M can be regarded as a measure on K via the projection 
K K/M: for / G C{K), 



^o(/) = / / f{km)dmdvQ{k) 
J K/M Jm 

for the probability invariant measure dm on M. Using this, if dg = y^^^^dydndk 
is a Haar measure on G = NAK, for £ Cc{G/M), 

(25) /x(V')= / / ^p{g)y^y-''-^dydnduo. 

Jk Jan 

This measure is left T-invariant and descends to the measure on {T\X) = 
V\G/M. 
Then 

Theorem 9.3. For ip G Cc{T^{T\X)), 

'4>{nay)dn ~ c^^jl{ip)y^~^ 



n&{NMr\T)\NM/M 



as y ^ 0. 
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Proof. See |23] . The argument given there works for general rank one space. 

□ 

10. Orbital counting on a light cone 

In this section we give a quantitative estimate of the orbital counting on 
]pn+i using the results in previous sections. 

Let r C G be a geometrically finite torsion- free subgroup with 6 > D/2. 
Let {vqjVq) = where ( , ) is the standard Hermitian form of signature (n, 1) 
on F"'-*^. Let L be the light cone of F"'^ given by {v,v) = 0. Then G acts 
transitively on L and so there is go £ G such that the stabilizer of uo5o is 
NM. By conjugating T hy go, we may assume that the stabilizer of vq in G 
is NM. 

In the language of section HI then vq is proportional to (0, 0, 1) and the 
action of Oy on it is 



1 



y 
10 
l/y 



I 



I.e., 

1 

VoOy = -Vq. 

y 

Traditionally we prefer G to act on the right and F on the left. So vectors 
in F"'^ will be row vectors and the matrices act on the right. In right action 
notation, vqNM = vq, so the stabilizer of the row vector vq is NM, and 
voay = ^vq. 

Lemma 10.1. If vqV is discrete, then the image of the horosphere Hq corre- 
sponding to Vq is closed in T\H^. Furthermore the image of the horosphere 
Hq is closed in T\Hf if and only if Hq is based either at a bounded parabolic 
fixed point or at the point which is not in the limit set. 

Proof. Any horosphere is obtained by projectivizing 11 intersection with the 
hyperboloid where 11 is an affine hyperplane Hx = {Y G F"'"'"^ : {{Z,X)) = 
1} with X G L\0. Here (( , )) is a standard positive definite inner product 
of signature n + 1 on F""*"^ and X is uniquely determined up to an element in 
C/(l) or Sp{l). This gives a homeomorphism between the set of horospheres 
and the set of {L \ 0)/U{l) {resp. Sp{l)). 

Now suppose the orbit vqT is discrete, and Hq to be the horosphere cor- 
responding to Vq. If T\Hq is not closed in T\H^, there are points pi G Hq 
so that Ti{pi) accumulates to a point p ^ T\Hq where vr : H^ — )• T\Hf. This 
means that there exist 7^ £ L so that piji accumulates to p G H^. Hence 
IfoTi accumulate, which is forbidden by discreteness of HqT. The second 
claim is proved by Dal'Bo [H] for the Hadamard manifold whose sectional 
curvature is bounded above by —1. She showed that the projection of the 
horosphere based at 00 is closed if and only if either 00 is a bounded par- 
abolic fixed point or 00 does not belong to the limit set provided that the 
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length spectrum of T is non-discrete, which is shown in [22\ for rank one 
symmetric space. □ 

Define a function on r\G by 

Prig) = ^ XBrivoig), 
■ye{rnNM)\r 

where Bt = {v £ vqG : \\v\\ < T}. Specially Fr(e) is the orbital counting 
function of T. Note that for ip G Cc{T\G) 



7e(rnAfM)\r 

/ [ X] XBt(W5)^(75)]'^S' is r invariant) 

XBT{vog)tp{g)dg ■ 



'{rnNM)\G 

Write G as G = {NM/M)AM{M\K) and accordingly g = [n]aymk. Then 
the above integral is 



Jm\K J\\voayk\\<T J[n^]£{rnNM)\NM/M 
'I\K Jy>T-^\vok\\ J [n^]e(rnN M)\N M/M JmeM 
J M\K Jv>T-'^\\vnk\\ J \n^]G(VnNM)\NM IM 



'M\K Jy>T-'^\\vok\\ J [n^]€{rnNM)\NM/M 

(26) = / / i;^{ay)y-''-'dydk, 

Jm\K Jy>T-'^\\vQk\\ 

where tpk{g) = J^^m ip{gmk)dm and ip^ {ay) = J(^rnNM)\NM/M i'i[n]ay)d[n] 
as before. 

Lemma 10.2. For any tp £ Cc{T\G), as T —)■ oo 

(Ft, V) ~ 5-'c^,T' [ fii^k)\\vok\\-'dk. 

JM\K 

If ijj £ C^{T\G)^ and || • || is K -invariant, 

{Ft,^^) = {i^,^o)S-'c^,T'\\vo\\-'{l + OiT-'')). 

Here the implied constants depend only on Sobolev norm of tp and and 
supp{ip). 
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Proof. The function tpk is M-invariant, hence it is defined on r^(r\X), and 
by Theorem 19.31 



ipk (ay) = / i;k{[n]ay)d[n] ~ c<i,^,fi{iJk)y 

Jne{rnNM)\NM/M 



as y — )• 0. By the equation (j26|) 

Jm\K J y>T-^\vok\\ 

= S-'c^,T^ [ fi{iPk)\\vok\\-^dk. 
Jm\k 

If ip and II • II are JC-invariant, since ipk = by Theorem 19.2^ 

^p''{ay) = {^^J,cPo)c^,y''-'{l + 0{y'')). 

Hence 

(Ft,V)= / / {^P,cPo)c^,y''-^l + 0{y''))y-''-'dydk 

Jm\K Jy>T-^\\vok\\ 

= dk{M\K) [ {^P,cPo)c^oy''~'i^ + Oiy''))y-''-'dy 

Jy>T-mvo\\ 

= dk{M\K){i,,cl)o)c^,6-^T%o\\~^{l + 0{T-^')). 

□ 

Theorem 10.3. As T ^ oo, 

Frie) ~ c^,6-^T^ [ 
Jk 

If \ \ ■ \ \ is K -invariant, 

Erie) = rVo(e)c<^olboir^T^(l + 0{T-^')). 

Proof. For smah e, choose a symmetric e-neighborhood of e in G, which 
injects to r\G, such that for all T » 1 and < e ^ 1, 

BrUe C and 5(i_,)t C r\u(iU,BTU. 

Let (/>e G C^{G) be a nonnegative function supported on Ue with cf)^ = 1. 
Define <l>e on T\G by averaging over T 

By definition, ^e{[g]) + only if \g\ G r\rf7e. Also for g G [7^, if ||z;o7|| < T 
then ||7;o75ll < + 
So 

^^'^ 7e(rnAfM)\r '^^ 7e(rnAfM)\r 



IwofeN-^dA:. 
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ye{rnNM)\r'^^^ 

Note that /^^^ XB(i+,)T(^o75)0e(5)c^5 < Iu,^<^i9)dg = 1 and the equahty 

holds only if \\vQ^Supp{(j)f,)\\ < (l + e)T. But if ||fo7|| < T, then ||fo7f^e|| < 
(1 + e)T. Hence 

Similarly 



7e(rniVA/)\r ' 

Since \\vQ^Supp{>P,)\\ < (1 - e)T if ||7;o7|| < T, < Frie). 

So we obtain 

(27) (i^(i-e)T, ^e) < Frie) < {F^i+,)t, ^.)- 

By Lemma 110.21 

~<5-%o(ni±e))' / iim)k)\\vok\\-'dk. 

J M\K 

But 

K{^e)k) = I I ^e{9mk)dmdfi{g). 
Jr\G Jm 

Hence 

\vok\\-^dk 



{F^,^,)T,'^e)^S-'c^,{T{l±e)Y [ [ [ ^,{gmk)dmdfi{g)\\, 

Jm\k Jger\G Jm 

= 5-^c^,{T{l±e)Y f [ <^,{gk)\\vok\\-'dkdil{g). 
Jr\G JK 



For a fixed k^, using Equation (I25p 

/ ^,{gko)d(i{g) = f ^,{gko)y'%dvo = [ M9)y'^iko'rduo. 
Jr\G Jr\G dk Jr\G dk 

But if e is small, y runs, say from 1 — e to 1 + e. Since J 4>edg = 1, 
Ir\G'^^(9)y^^{ko^)*diyo is in the order of f^^^y^ = 0{e). 
Hence the above equality becomes 

<!>,) ~ ric^o(T(l ± e))'{ [ \\vok\\-'dk + 0(e)). 

JK 

In equation ([27|) . let e — )• to obtain 

\vok\\-'. 



JK 

nt, take and ET- 



IK 

If II • II is ET-invariant, take Uf and (h^ ET-invariant, and we obtain 



EQUIDISTRIBUTION 



31 



Since must contain K and dm is a probability measure on M G K, we can 
take (peie) = 1 to have / (p^dg = 1. Hence {^e,(po) = Jg(zu,^49)(l>oi9)dg = 
00 (e) + Re and i?^ — as e — )• 0. So we obtain 

Frie) = Me)c^o^-HT)'\\vo\\-'il + OiT-'')). 

□ 

Finally the orbital counting theorem can be stated as follows: 

Theorem 10.4. Let T C G = KAN be a geometrically finite group with 
critical exponent <5 > y for real, complex and quaternionic hyperbolic space, 
where KAN is a fixed Iwasawa decomposition introduced in sectionlM Sup- 
pose vq is in the light cone such that vqT is discrete, and the stabilizer of vq 
in gQ^TgQ is in NM. Then for any norm || • || on F"'""^ 

#{v G ^;or : \\v\\ < T} ~ c^^S'^T^ [ \\vo{go^kgo)\\-^dk. 

Jk 

//II • II is gQ^ KgQ-invariant, then 

#{v G voT : \\v\\ <T] = c^,5-^T%o\\-\l + 0{T-^')). 

11. APPLICATION TO APOLLONIAN PACKING 

In euclidean n-space, the maximum number of mutually tangent n — \- 
spheres is n + 2, if one demands that no three spheres have a common point. 
An Apollonian sphere packing of n-dimensional euclidean space is for a given 
n + 1 mutually tangent spheres, one adds a sphere which is tangent to all 
the n + 1 previous spheres. This procedure defines a sphere packing, called 
Apollonian sphere packing. 

The formula for curvatures for such an n + 2 mutually tangent spheres is 

mm 

n+2 n+2 
1=1 i=l 

For a given n + 1 mutually tangent spheres, there are exactly two spheres 
satisfying the above equation. The sum of curvatures of these two last 
spheres is 

2 

n — 1 ^-^ 

i=l 

Hence for n < 3, if the first n+2 generating spheres have integral curvatures, 
then the curvature the next sphere is also integral. So we restrict 

our attention to integral Apollonian packings. 
Let 

n+2 n+2 
Q{ki,--- ,K„+2) =n'^Kj - Kj 
1=1 i=l 
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be a quadratic form on M""*"^. This has signature (n + 1,1). Hence the 
orthogonal group Oq can be identified with the group of hyperbolic isometry 
in H^+\0{n+ 1,1). 

Prom now on we focus on = 3, hence 

4 

1=1 

There is an integral group, called Apollonian group, defined as follows. 

Changing the curvature - ,K4,K5) to (ki,--- ,K4^,k'q) corresponds to 
the integral matrix 

/l \ 

10 

^5= 10 

1 

Vi 1 1 1 -V 

Similarly there exist corresponding matrices Si, S2, S3 and ^4. Denote A = 
{Si, S2, S3, S4, S5) C Oq{Z), the Apollonian group. This group acts on 
apoUonius sphere packing V generated by the initial spheres with curvature 
Ki, - ■ ■ , K5 . We will see shortly that this Kleinian group is geometrically 
finite. 

If S^, - ■ ■ ,S^ are initial spheres, then next sphere S^ is obtained by a 
reflection along the sphere which passes through the intersection points 
of S^ , S'^ , S^ , S'^ . So the element ^5 corresponds to a reflection along T^. 
It is similar for other elements Si. This way A corresponds to a group F 
generated by reflections 

Let Di be a hemisphere in upper half space of R'^, whose boundary is 
T*. Since is the ideal boundary of and the reflections in gener- 
ates whole isometry group of H^, we know that F is a discrete subgroup 
of Iso{H^) generated by reflections along Di with a fundamental domain 
bounded by -D/. Since the fundamental domain has finite sides, it is a geo- 
metrically finite group with cusps. 

Note that F acts on V by permuting spheres in V, so it leaves invariant the 
set of tangent points between spheres. But it is easy to see that any point 
on some sphere can be approximated by tangent points, hence the limit set 
of F is equal to the union of spheres in V. Since a sphere has Hausdorff 
dimension 2, the Hausdorff dimension of the limit set Ap > 2 > 3/2. Hence 
our assumption that the critical exponent is greater than D/2 is satisfied to 
apply previous results to this section. 

11.1. Orbital counting for A. We use a ivT-invariant norm || || in this 
section, in fact a maximum norm, and G = S0'^{4, 1). As before Bt = {v e 
: Q{v) = 0, ||f|| < T}, i.e., the set of curvatures of maximal mutually 
tangent spheres in whose norm is less than T. Fix an initial curvatures 
^ so that (5(0 = ^-iid the corresponding spheres generate an Apollonian 
sphere packing in R^. By conjugation if necessary, we will assume that the 
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stabilizer of ^ is equal to NM. One has to recall the definitions and notations 
from section [lOl Let c/.^ S C^{H^) = C^{SO^{A,l))^ be a nonnegative 
function supported on a if- invariant with Jj^4 (p^ = 1. Define a function 

defined on M = r\H^ = T\G / K by 

Fix a primitive integral polynomial / in 5 variables. One defines a positive 
sequence A{T) = {an{T)} where 

aniT) = WT{-i) 
7eStabr(C)\r,/(e7)=™ 

where wxij) = Jq/j^ XBT{^19)4'ei9)dfJ'{g) for each 7 G F. A subsequence 
Ad{T) = {an{T)\n = O(mod d)} and define \A{T)\ = ZMT), \Ad{T)\ = 

Z]n=0(mod d) (^n{T). 

For any subgroup Fq C F with Stabr(C) = Stabro(C), define 
4° (5)= E XBriCig), 

7eStabr{0\ro 

and 

<7i(5) = E "^^(^r'^^), Tier, 

7ero 

which is an e-approximation to the identity around [7^"^] in To\H^. Finally 
let 

F5(d) = {7 E F : e7 = e(mod d)}, 
and its finite index subgroup 

T{d) = {7 G F : 7 = /(mod d)}. 

Then Stabr(C) = Stabr^(rf)(C). The following is due to Bourgain-Gamburd- 
Sarnak [5], Salehi Golsefidy-Varjii [32j, Breuillard-Green-Tao [7J and Pyber- 
Szabo [3T| . 

Theorem 11.1. For any square free integer d, the -spectrum ofT{d)\H^ 
has a uniform spectral gap 2 < 6 < 5r so that it does not have any spectrum 
in [e{3-e),6ri3- 5r)). 

Proof. (Step I; Family of expanders) For general Apollonian sphere pack- 
ing in M", 

r = A={Si,--- ,Sn+2) cG = Oq{Z[^—]) C GL{n + 2,Z[^—]). 

n — 1 n — 1 

Note that F is not discrete if n > 4. Then the following is proved by A. 
Salehi Golsefidy and P. Varjii [32] , 
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Theorem A Let T C GL{d, be a group generated by a symmetric 

set S. Then CayleyiV /T{q); S mod q) form a family of expanders when q 
ranges over square-free integers coprime to qo if and only if the connected 
component of the Zariski-closure of T is perfect. 

(Step II; Uniform spectral gap) This is a corollary of the above theo- 
rem. For reader's convenience, following closely the argument in [5], we give 
the details of its extension to higher dimensional real hyperbolic spaces. 
Suppose r is a geometrically finite discrete group in G = 0(n, 1). For large 
enough q, by [26] T/T{q) = G{¥q) where Fp = Z/pZ. 

Fix a fundamental domain T C H^'^^ of F. For any map / defined on 
T{q)\H^+\ let / be a lift to H^+\ Then / can be regarded as a vector 
valued function F defined on T\H^'^^ by setting 

I^iz) = (/(72))7er/r(g) 

where z G by identifying T\H^~^^ with T. It satisfies the equivariant 
condition 

F{jz) = R{^)F{z) 

for 7 G G(Fg), where R{'^) denotes the right regular representation of G(¥q). 

More formally, every L? function on T{q)\H^^ can be viewed as an L? 
section of a flat bundle E of rank |G(Fq)| on r\i?j^'*"^ as follows: Given 
/ G L\T{q)\Hl+\^), define F : H^+' ^ L\G{¥q)) by 

F{x,g) = f{gx) 

for X G g G F. Then, for /i G F, 

F{hx,g) = f{ghx) = {R{h modF(g))F)(x, 5), 

so F can be viewed as a section of the flat vector bundle E over T\H^^ 
associated to the F-principal bundle H"^"^ — t- T\H^^ via the right regular 
representation of G(Fg) composed with F — )■ F/F(g) = G{¥q). Or, equiva- 
lently, E is associated to the G(Fg)-principal bundle T{q)\H^^^ r\i?^+^ 
via the right regular representation R of G{¥q). 

Then E = Eq + Ei according to the splitting of the right regular represen- 
tation into constant functions, and functions on G{¥q) which are orthogonal 
to constant functions. This splitting is orthogonal and parallel (compatible 
with the flat connection). The splitting of F = Fq + Fi corresponds to 
/ = /o + /i where /o is G(Fq)-invariant. 

If / is orthogonal to the Aq eigenspace of T{q)\H^'^^, which is generated 
by a G(Fq)-invariant function, then each factor of / is orthogonal to the 
Ao eigenspace of F\ifj^~'"^. Integrating each factor over T\H^'^^ and taking 
summation, we get HV/lH > Ai(F\i/^+^)||/|||. Here || • II2 denotes the 
L^-norm over T{q)\H^'^^. 

Hence from now on, we can assume that f = fi, F = Fi takes values in 
El. 
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Pick a generating system S of F. Construct corresponding right Cayley 
graph. By definition of the spectral gap for the discrete Laplacian, for each 

^ Y^mz) - R{s)F{z)\f > Ai(G(F,))||F(z)||2. 
Integrate over r\ff]^~'""^ and get 

l"^! ses 

Let Hq be the subspace of such functions orthogonal to the bottom eigen- 
function cpo. We need to show that there exists e > independent of q such 
that for any F in Hq 



I^\\VF\\^dfi 



> An + e. 



Above discussion about expanders implies that for z £ J-, and for any 
F £ Hq, there exists 7 G such that 

(28) ||F(7z) - F{z)\\' > Ai(G(F,))||F(z)||2, or 

(29) ||i?(7)F-F||i> Ai(G(F,))||F||2. 

Set / = I |F|| = a(l)Q{z)+b{z) where b is orthogonal to (j)Q, the Aq eigenfunction 
normalized that Jjr<pQ = 1, and norm of F, + J-^b{z)^dfi is 1. Then 



one can show that 



J^\\VF\\^dfi 



> Ao + (Ai - Ao) / b^dfi. 



If fjr b'^dfi is bounded below for all F G Hq we are done. Hence suppose 
there is no lower bound. Then after taking a weak limit, we may assume 
that there is F G Hq with J-^b'^dfi = and a = 1. One can write 

F = {Fi), ^ = Ui 

so that "^uf = 1, which implies that "^Uj^^ = 0. Set u = (ui). Then 
Equation (j28|) becomes 

(30) ||i?(7)^(z)-n(z)|p>Ai(G(F,)) 

for any z £ T and some 7 G 5. By Step I, Ai(G(Fq)) is uniformly bounded 
below independent of q. A direct calculation as in [S] shows that 



|V(0oti)ll'^ = X3 |V(0|,U;i)|- 
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which imphes that 

J^\\VF\\^dfi _ /^||V(/>o|p + 0g||V^|pd// ^ ^ J^(t>l\\Vu\\^dfi 

If is bounded below, we are done again. Remember that we 

normahze (po so that Jjr4'Q = 1- 

Now u = (ui, • • • , Uk) for k = \G{¥q)\ and hence Vn = (Vui, • • • , Vuk)- 
Fix zo G Take a unit speed geodesic a{t) connecting zq and 720. Then 

u{jzo) - u{zo) = I ^M^^(it= I Vu{a{t))a{t)dt. 



Hence 



|n(7^o)-^(^o)||' = (/ ^^^^^^dtf < d{zo,7Zo) I \\Vu{a{t))-a'it)\\'dt 



<d(zo,7^o) / \\Vuiait))\\'dt 

J a 

Take a Fermi coordinate along a so that metric tensor is 

ds^ = dp^ + sinh^ pd(j)^ + cosh^ pdt^ 

where p is a distance from a so that the volume form dp = cosh p sinh"~^ pdtdpdcj). 
Take a small cross section B orthogonal to a passing through zq so that the 
induced volume form on B can be written as dB = sinh"~^ pdpd(f>. Then 
N = B X a will be a small tubular neighborhood around a so that 

L = max{d(2;, 72;) : z G B} < 00, m = imncpQ > 0. 

Since (po is a F-invariant positive function, m > 0, see section [8l By inte- 
grating above inequality over B, 

\u{-fz) -u{z)\\^dB < L \\Vu\\^dtdB<L \\Vu\f cosh pdtdB 

Jn Jn 

= L [ \\Vu\\'^dp < — [ cplWVuW'^dp < — [ (plWVuW'^dp. 
Jn fn Jn rn Jjr 

By equation (1301) . J -p (jy^WV u\\'^ dp has a uniform lower bound. □ 

Since F(d) is a finite index subgroup of r^(d), such a spectral gap theorem 
holds for T^{d) as well. Prom this and Lemma 110.21 we obtain 

Proposition 11.2. There exists eo = 5' uniform over all square free integer 
d so that for any 71 G F and for any V^{d), we have 



B 
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Proof. By Lemma[in21 for ^1^}^^ £ C^(T^{d)\G)^ and K-invariant || • ||, 
since T^{d) is finite index in F, 5r5((i) = <^r and 

Here since 6' depends only on the spectral gap and since it is uniform over 
all square free integers d, it is uniform. Let (f)o be a lift to M{d) = r^{d)\H^ 
of the bottom unit eigenfunction 0o on M = T\H^. Then the unit bottom 
eigenfunction on M{d) = T^{d)\H^ is 

r id) 

Since $e,7i is an e-approximation to the identity around [7""^] in M{d) and 
(j)Q is invariant under F, 



1 , > 1 



= W^f'^'^'^ '''''' 

Also note that c r,(d) = ^ , csn from Section [5l Here once we fix 6^, 

the implied constant depends only on Faj,^^^^. But since r^(d) is of finite 
index of F, it is constant. Hence the claim follows. □ 

Corollary 11.3. There exists eo uniform over all square free integer d such 
that 

01(d) 

where Oj{d) = J2^^^r^{d)\r, /(C7i)=o(d) I and x = 5^^c^,d,\\C\\-^^T^^ . 
Proof. 

n=0{d) 7eStabr(0\r 7ierj(d)\r 7GStabr(C)\r5(<i) 

/(«7)=0(rf) /(57i)=0(d) 

7ier5{d)\r, /(C7i)=oM 7estabrK)\re{rf) '^/^ 



5^ / F^«^')(5)0e(7r^ff)rf^(5) 



7ier^(d)\r, /(^7i)=o{d) 



7ier5(d)\r, /(57i)=o(d) •^r5{d)\G/i^ 
7ierj(d)\r, /(?7i)=o(<i) 
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By Proposition 111.2^ 

7iere(d)\r, /{f7i)=o(rf) ' 

- (^ + 0(r5r-eo))^ 



[r : T^{d)] 

where Oj{d) = E^^(,r^{d)\r, /«7i)=o(d) 1 and x = S^^c^^d,\\C\\-^^T^^ , hence 
the claim follows. □ 

11.2. Digression to an algebraic group. Let G be a semisimple alge- 
braic group defined over Q. Denote G{¥p) a reduction of G to an alge- 
braic group defined over ¥p = Z/pZ. In our case G = Spin{A^ 1), a dou- 
ble cover of S0^{4:,1). Note that a stabilizer of a point ^ ^ V = {v = 
(xi,--- ■,x^)\Q{v) = 0} \ is a connected Q-group H = NM, so that 

V = G/H. By [30] (Proposition 3.22), the reduction of V to V{¥p) = 
r5 



{v = {vi,--- ,V5) e : Q{v) = 0{p)} is G{¥p)/H{¥p). By [26], since T is 



■p 

Zariski dense, the reduction map 

r ^ G{¥p) 

is surjective. Note that T^{p) is a stabilizer of ^ in V(¥p) where is a 
p-reduced image of ^ in V{¥p). Since the reduction of F is G{¥p), 

by m- 

Let /i(xi,-- - ,X5) = xi, /2(a;i,-- - ,^5) = 3:1X2. To estimate Oj^^^^ = 

^7i6r5(p)\r, /i(57i)=o(p) 1, if Q'(2;2, • • • ,X5) = Q{0,X2, - ■ ■ ,X5) is a quadratic 
form whose zero set is = {(0, X2, • • • , X5) G y}, then 

Ol(,. = i^Wi¥p)^p' 



Mp) 

by [2]. Similarly 

"''f2{p) 

Then 



(^l(v) = *W{¥p) ~ 2/. 



In general, if fk{xi,--- ,2:5) = xiX2---Xfc, then gk{p) ~ A;p ^ for A; = 
,5. 

11.3. Asymptotic growth of number of spheres with prime curva- 
tures. Let /i(xi,--- ,X5) = Xi, /2(xi,-- - ,X5) = XiX2,/fc(xi, • • • ,X5) = 

X1X2 • • • Xfc as before and for square free integer d, let 



9i{d) 



Olid) 



[T:r^{d)r 

Proposition 11.4. There exists a finite set S of primes such that 
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(1) for any square free integer d = did2 with no prime factors in S and 
for each i = 1,2, gi{did2) = gi{di)gi{d2) 

(2) for any prime p outside S, gi{p) G (0,1) and gi{p) = p~^+0{p~^), g2{p) = 
2p~^ + 0{p~'^) for some q > 0. 

Proof. The second claim is already shown. By [26], for a large prime p, the 
reduction of T is G(¥p). So let S be the set of primes which are less than 
such a prime p. Then for d = pi ■ ■ ■ pk square free with pi ^ S, the diagonal 
reduction 

r ^ G{Z/dZ) G(FpJ X • • • X G(FpJ 

is surjective and it follows from Goursat's lemma that F surjects onto G(Z/diZ) x 
Gil^jd^I^ for any square free d = d\d2 without any prime factor in S. Hence 
for d = did2 square free without prime factors in S, the orbit of ^ mod d un- 
der r is equal to the orbit of G{Z/diZ) x ^(Z/daZ) in {Z/diZf x (Z/diZ)^ 
The same thing is true for the orbit satisfying the equation /(^7) = 0{di). 
Therefore g{d) = g{di)g{d2). See [6l|23]. □ 

Note that a„(r) = E7eStabr(«)\r, /i(g7)=n ^t(7) is a smooth counting 
function for the number of vectors (xi, • • • , 2:5) in the orbit ^A^ of maximum 
norm bounded above by T and xi = n. We want to obtain an asymptotic 
growth of the number 7r'^{T) of spheres in an Apollonian packing V whose 
curvatures are prime and less than T. The initial spheres with curvature 
^ = (^1, • • • ,.^5) is chosen so that ,^1 < which corresponds to the largest 
sphere bounding all the other spheres, and (,2,(,3,S,ijS,5 the smallest curva- 
tures (largest spheres) in the packing. Upon iteration, one obtains a smaller 
sphere, hence larger curvature. 

Since any sphere in the packing is obtained by the initial five spheres 
corresponding to ^ under the group A, with the maximum norm on M^, 

7r^(r) < 4 + #{7 G ^ : 1 1^7*1 Uax is prime < T} 

5 

< G : ll^llmax < T, Vi is prime}. 

i=l 

Similarly the number tt2{T) of twin prime curvatures less than T is 

TT^{T) <^ 7^(7 G A : ||C7*||max is prime < T, one more entry of ^7* is prime} 

5 

< ^^#{1; = (fi,--- ,^5) G (,A^ : ||w||max < T, Vi,Vj primes}. 

This counting problem uses so-called a Selberg's sieve; using S in Propo- 
sition 111.41 and g as one of (7,, a multiplicative function h on square free 
integers outside 5 that h{p) = i^^g^p^ for a prime p outside S, 

Theorem 11.5. ([20], Theorem 6.4) Let P be a finite product of distinct 
primes outside S such that for any square free d\P 

\A^{T)\=g{d)x + rdiA{T)). 
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Then for any D > 1, 

S{A{T),P)= Y «n(T)<x( E Kd))-'+ Yl Md)\r,{A{T))l 

(n,P)=l d<VD,d\P d<D,d\P 

where r^id) denotes the number of representations of d as the product of 
three natural numbers. 

Prom all these we obtain 

Theorem 11.6. Given a bounded primitive integral Apollonian sphere pack- 
ing V inM.^, 

vr^m « ^ 
logT 

and 

Generally ifTT^{T) denotes the number of k-mutually tangent spheres whose 
curvatures are prime numbers less than T , then 

for k < 5. 
Proof. Since 

\A,{T)\ = .f'^'^l.A x + 0(r^r-°)) = gi{d)x + rd{A{T)), 



and since < 1, rd{A{T)) < T-^r-eo. For any ei, EdKoMd) < 



r:T^{d) 

1^ + 2^ H h -D^ << D^+'S hence 

Y ■r:i{d)\rd{A{T))\ < D^+'^T^^-'° < T^^/logT 

d<D,d\P 

if D = j"=o/4_ jjence take P to be the product of all primes less than T*^"/^ 
outside S. Denote = 1 if n is square-free and otherwise. Also let 

s{m) denote the largest square-free number dividing m. Then for h = jr^; 
one can deduce that ([20], Section 6.6) 

Y h{d)r. Y n^MT^ 

d<VD,d\P d<y/D,d\P 

= E ^~'^p\^T^= ^ ^"'n,|,(i + i + l + ...) 

d<y/D,d\P d<^/D,d\P 
V- .111 ^/. 1 1 N 1 1 

^ d dpi dpi P2 Pi Pk Pi 

d=pi---pk<VD,d\P 
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But since the numbers appearing in the denominator is of the form m = 
Pi^ • • • and since 5 is a fixed finite set of primes less than some number, 
and P is the product of distinct primes less than D outside S, when D tends 
to infinity, the above sum is 

» E (Mrf) E -) = E - 

d<^/D m,s{m)=d m,s{m)<VD 

> J2 > log D > log T. 

m<\/D 

Hence Theorem If 1.51 gives 

rpSr 

siAinp)«^. 

Since there are only 5 coordinates, we may assume that for any v = 
{vi,--- ,^5) G ^A^ vi is the largest coordinate among Vi. Therefore by 
Equation (j27p for the product P of distinct primes less than outside 
S, 

-— »5(^((f + e)r),p) + (r^o/Y^= V a„((f + 6)r) + (r^o/Y^ 

^ (n,P)=l 

» #{V = (xi, • • • ,X5) G -Ml max < T, [x^, P) = 1} + [T'^/^f^ 

» (r.o/4)5r + ^{^ = (3.^^ . . . ^^^) e ^^4* . 1 1 < T^^^ prime > 

> 5#{i; = (xi, • • • ,X5) G : \\v\\ max < T,xi prime} 
> = (^^i,--- ,2:5) G : ||w|| max < T,some Vi prime} > vr^(r). 

Let u;{d) denote the number of distinct prime factors of d and T^ld) the 
number of positive divisors of d including 1. Then T^id) = 2'^^^^ for square 
free integer d. For g2, considering is a finite set and d\P implies that d is 
a square free, as before 

E ^('^)~ E "pM i!^2p-i 

d<vl),d\P d<'/D,d\P 
d<VDA\P d<VD 

= y: ^^{d)v{d) y: i- 

d<\/D m,s{m)=d 

If m = ■ ■ ■ p'j^'' for distinct primes pi then s{m) = d = pi- ■ -pk- Then 
V{d) = 2^ V{m) = (ai + 1) . . . [at + 1). Hence 

V{s{m)) T>{m) 
d m 
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For each m < y/D, if s{m) = d < m, then appears in 

d<^/D m,s{m)=d 

Hence 

P(m) 



E Mrf)i^(rf) E > E 

Note by [27] (page 94), 



m ^ — ' m 

d<v^ m,s(m)=d m<VD 



J2^ = l(^ogNf + OilogN). 
^-^ n 2 

n=l 

From this 

M^^)»aogZ))'»(iogr)2. 

d<^/D,d\P 

-^^»5(yl((l + 6)r),P)= Yl (^ni{l + e)T) 

(n,P)=l 

> #{v = (xi, • • • , xs) G : max < T, {xiX2, P) = l} + (r^o/4)^r 

» #{V = (Xi, • • • ,X5) G : Ibll max < T,Xi,X2 prime > r^o/4| + (^eo/4^5r 

»vrr(T). 

For general A; < 5, it suffices to estimate 



E ^PMr 



d<vl5,d\P d<\/D 



= Kd)k-(^^ E -• 

d<VD m,s{m)=d 

If m = p"^ • • • for distinct primes pi then s(m) = d = pi ■ ■ ■ pk- Hence 

d ^ m 

For each m < \/^, if s{rii) = d < m, then ^"^''^ appears in 



Hence 



m 

d<VD m,s(m)=d 



E ^(d)k-''' E ;!;> E 



m — ' m 

d<VD m,s{m)=d m<y/D 
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where the last inequaUty follows from the following theorem [36], and com- 
municated to us by Pieter Moree. 

Theorem 11.7. 

^ fc^W = cfc(0)x(logx)'=-i + 0(x(logx)'=-2), 

n<.x 

where 

□ 

Indeed, using a finer version (Theorem 11.13 of [13J) of Sieve method of 
a multiplicative function g{p) = | + 0{p~^~^) 

Theorem 11.8. There exists a real number I3{k) > such that 

if is) + 0(log D)-i/6)|^(T)|np|p(,)(l - g(j,)) + R{D) < 5(^(T), P(z)) 

< (Fis) + 0((log D)-V6)|^(r)|np|p(,)(l - g{p)) + R{D) 

where z = D^/^ with s > ji{k), and P{z) is the product of distinct primes 
less than z outside S, where F{s) > and /(s) > are certain functions of 
s > 0, depending on k, defined as solutions of explicit differential-difference 
equations, such that 

lim /(s) = lim F{s) = 1, 

and where 

R{D) = Y,\rM{T))\. 

d<D 

In both upper and lower bounds, the implied constant depends only on k and 
on the constants in the asymptotic of g{p). 

we can try to give a lower bound for our counting problem. Note that the 
number vr^ (T) of spheres whose curvature is prime less than T is equal to 

#{u = (xi, • • • ,X5) e : \\v\\ max = prime < T}. 

But an{T) approximately counts the number of orbit vectors with maximum 
norm less than T and whose first coordinate is n. Hence as long as the 
first coordinate is prime p, both {p, V2,V3,Vi,V5) and {p, V2,v'^,v'^,v'^) will be 
counted even though none of ^2,^3,^4, is prime. This means that OniT) 
with n prime, will count the same sphere as many times as it appears in the 
orbit ^A^. Another difficulty is that 

#{^; = (xi, • • • , X5) G : I I max < T, T^^/^ < xi prime < T} 

is not equivalent to 

#{t; = (Xi,--- ,X5) e : llvll max < (Xi,n . ^^lop) = l}. 

nrime v<Ti 
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To make them comparable, we pose the condition that xi cannot be written 
as a product of more than r primes > T^o/^ where r is chosen as the first 
positive integer to satisfy T 4 > T. People call it r-almost prime. 

These things complicate the problem for the lower bound. But if we allow 
this over-counting with extra assumptions, one can give a lower bound. 

Theorem 11.9. Let 7r^{TY denote the number of 5 spheres kissing each 
other (i.e. the number of orbits) among whose at least k curvatures are r- 

almost primes and all of whose curvatures are less than T where r is the 
first positive integer larger than ^ for a fixed large sq. For any k <b, 

Proof. We know already that |^(r)| T^^ . Note also that since rd(^(r)) < 

R{D) = Y,\rd{A{T))\<DT'^-'\ 
d<D 

Hence if D = T^o/^, R{D) < (i^^- Now we estimate np|p(^)(l - g{p)) 
using Mertens formula, which is communicated to us by Pieter Moree. Put 
Akip) = l — k/p\ip> k and A}.{p) = 1 otherwise. Put = Wp>2 Ak{p){^ — 
\/p)~^. Let 7 denote Euler's constant. 

We claim that as x tends to infinity we have Wp^xAkiv) ~ Bke~^'^ \og~^ x . 
Write 

n Ak{p) = 11(1 - i/pf n Ak{p){i - i/pr'. 

p<x P<x P<x 

Note that log^fc(p) - A:log(l - l/p) = 0{k'^/p^) as p tends to infinity. It 

thus follows that the latter product converges to as x tends to infinity. 
Noting that 

I lop- /li.lrjl — A- Inp-I I — I /nil = f ;( > ' ' 



I \ogAk{p) - /clog(l - l/p)\ = 0{Y^ ^) = O(-), 



X 

p>x p>x 

wc find that 

n Ak{p) = - vpfBk{i + o(^)). 

p<x P<x 

Using the Mertens theorem that 

n(i-vrt~i^. 

p<x 

the claim then follows. 

Now take D = T^o/i ^nd let s large enough and e ^ to get for z = Tia 

eg 

and the product P = P{z) of distinct primes less than outside S 

rpSr 

(a)HlogT)fe « {f{s) + 0{logD)-'/')\A{T)\Up\p^,){l-9{p)) + R{D) 
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«5(^((l-e)r),P) 

= XI «n((l-e)r) < = (Xi,--- ,X5) G : max < r, (xiX2---Xfc,P) = 1} 

(n,P)=l 

Since there are only 5 coordinates, there are only 5Ck choices of k coordi- 
nates out of 5, 

< #{V = {Xi,--- ,X5) e : max < 

5Ck#{v = (xi,--- ,X5) G : ||?;|| max < T, (xiX2 • • • x^, P) = l,31<i<k,Xi> t£} 

< = (xi,--- ,X5) G : ||w|| max < T,(xiX2---Xfe,P) = 1, Xi > TS,i = 1,... ,A;} 

+(rS)'5r 

^ #{'1' = (a^^i) • • • jX^) E : \\v\ \ max <T,T > xi,X2, • • • ,Xjfc r-almost prime > T*^} 

□ 

Note that to get an upper bound, we need to reverse the inequality 

#{V = (Xi, • • • ,X5) G CA^ : max < T, (xiX2 • • • Xfe,P) = 1, Xi > TS}+(TS)5r 

< #{v = {xi, ■■■ ,X5) G ^A : ||v|| max < T,T > xi,X2, • • • , Xjfc r-almost prime > T^,}^ 
which is not obvious to us. 

Corollary 11.10. Given a bounded primitive integral Apollonian sphere 
packing V in , Tr^{Ty denote the number of 5 spheres kissing each other 
whose curvatures are r-almost primes less than T where r is a fixed positive 
integer depending only on Apollonian packing. Then 
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